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ABSTRACT 

Measurements of the intensity of the cosmic X-ray background (XRB) carried out 
over small solid angles are subject to spatial variations due to the discrete nature 
of the XRB. This cosmic variance can account for the dispersion of XRB intensity 
values found within the ASCA, BeppoSAX and ROSAT missions separately. How- 
ever there are differences among the values obtained in the different missions which 
are not due to spatial fluctuations but, more likely, to systematic cross-calibration 
errors. Prompted by recent work which shows that ROSAT PSPC has calibration dif- 
ferences with all the other missions, we compute a bayesian estimate for the XRB 
intensity at 1 keV of lO.Oto'g keV cm ~ 2 s 1 sr 1 keV -1 (90 per cent confidence errors) 
using the ASCA and BeppoSAX data points. However, this value is still significantly 
larger than the HE AO — 1 intensity measured over many thousands of square degrees 
(8keVcm- 2 s- 1 sr" 1 keV -1 ). 
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1 INTRODUCTION 



Most of the X-ray background (XRB) above photon ener- 
gies of a few keV is known to be extragalactic in origin (see, 
e.g, reviews by Boldt 1987 and Fabian & Barcons 1992). Its 
spectrum was well measured by the HE AO — 1 A2 exper- 
iment in the 3-50 keV energy range (Marshall et al 1980), 
where it fits a thermal bremsstrahlung model with a tem- 
perature kT ps 40keV. An overall fit to the XRB spectrum 
from 3 keV to 10 MeV was presented by Gruber (1992) 
based on data from the A2 and A4 experiments on board 
HE AO — 1, which showed an extrapolated intensity at 1 
keV of » 8keVcm- 2 s" 1 sr" 1 keV -1 . 

More recently, measurements of the XRB spectrum 
and intensity at photon energies < 10 keV have been ob- 
tained with the use of imaging instruments. Using ASCA 
data, Gendreau et al (1995) confirmed that the thermal 
bremsstrahlung shape (which at photon energies below ~ 
20 keV can be approximated by a power law with energy 
spectral index a ~ 0.37) provides a good fit to the XRB 
spectrum down to about 1 keV. Vecchi et al (1999) used the 
LECS and MECS instruments on board BeppoSAX in the 1- 
8 keV band to confirm that the XRB spectrum is consistent 
with a power law with a ~ 0.4. 

ROSAT PSPC observations have also provided data 
on the spectrum of the extragalactic XRB. Hasinger (1992) 
found a very steep spectrum for the XRB below 2 keV 
(a ~ 1). However, shadowing experiments with ROSAT 
have provided the more stringent upper limits to the slope of 



the extragalactic XRB at soft X-ray energies (a < 0.7, Bar- 
ber & Warwick 1994). The relative steepness of the ROSAT 
spectra with respect to many other instruments has been 
studied in detail by Iwasawa, Fabian & Nandra (1999). The 
fact that the spectral shape of the same sources is usually 
consistent in observations of many instruments with the ex- 
ception of ROSAT (which usually finds steeper spectra) is 
highly suggestive of a calibration mismatch which will cer- 
tainly affect the XRB spectrum as well. 

The normalisation of the extragalactic XRB (that we 
parametrise as the XRB intensity Ixrb at 1 keV in units of 
keVcm -2 s _1 sr" ^keV" 1 ) remains uncertain. The Marshall 
et al (1980) measurement is the most robust result as it was 
performed over a very large solid angle (~ 10 4 deg 2 ) with 
instrumentation especially designed to subtract efficiently 
and accurately the particle background. The XRB intensity 
measured by the imaging instruments on board ASCA, Bep- 
poSAX and R OS A T invariably yield higher values which are 
often statistically discrepant among them. 

In this paper we point out that the discrete nature of 
the XRB introduces a cosmic variance in its intensity which 
is large for observations carried out over small solid angles. 
This cosmic variance is just the confusion noise caused by 
unresolved or non removed sources in the images (Scheuer 
1974, Barcons 1992). We find variations of the order of 10 
per cent for solid angles under ldeg 2 and therefore cosmic 
variance often dominates over the statistical uncertainties 
quoted in the various estimates of the XRB intensity. Once 
this is included, we combine individual measurements of the 
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Instrument Field Qi Ii Ui Sj Ref 

Name (dcg 2 ) (kcV cm~ 2 s _1 sr _1 kcV -1 ) 



ROSAT PSPC 


QSF3 





.223 


11.4 


0.34 


1.00 


Chen et al (1996) 




ASCA SIS 


QSF3 





.134 


10.0 


0.37 


0.76 


Chen et al (1996) 




ROSAT PSPC 


GSGP4 





.283 


11.8 


0.44 


4.49 


Georgantopoulos et al 


(1996) 


RObAl PSPC 


SGP2 





283 


12.0 


0.61 


4.53 


Gcorgantopoulos et al 


(1996) 


ROSAT PSPC 


SGP3 





.283 


12.0 


0.68 


4.52 


Georgantopoulos et al 


(1996) 


ROSAT PSPC 


QSF1 





.283 


9.9 


0.65 


4.49 


Georgantopoulos et al 


(1996) 


ASCA CIS 


Lockman 





.165 


10.9 


0.61 


1.31 


Miyaji et al (1998) 




ASCA CIS 


Lynx 





144 


9.3 


0.61 


1.38 


Miyaji et al (1998) 




ROSAT PSPC 


Lockman 





.165 


10.0 


0.30 


2.10 


Miyaji et al (1998) 




ROSAT PSPC 


Lynx 





.144 


11.5 


0.43 


2.28 


Miyaji et al (1998) 




ASCA SIS 


Various 





538 


8.9 


0.50 


0.96 


Gendreau et al (1995) 




ROSAT PSPC 


Various 


5 


9 


13.4 


0.18 


0.66 


Hasinger (1992) 




SAX LECS+MECS 


Various 


0. 


.726 


11.0 


0.30 


0.77 


Vecchi et al (1999) 





Table 1. Measurements of the XRB intensity at 1 keV carried out with imaging instruments. Qi is the solid angle covered, Ii the 
measured intensity, <t; the 1-sigma statistical uncertainty and the 1-sigma uncertainty derived when statistical and cosmic variances 
arc included. 



XRB to compute bayesian estimates of the overall XRB in- 
tensity. We find that measurements carried out within the 
same mission are brought to consistency by cosmic variance, 
but that systematic differences between different missions 
remain. 



2 THE COSMIC VARIANCE 

Table 1 lists the XRB intensity measurements that we have 
used. In particular and for each particular datum we list 
the measured intensity at 1 keV Ii, the quoted statistical 
1-sigma error Oi (converted from 90 per cent errors) and the 
solid angle covered Qi. It is clear that the variations between 
different measurements carried out within the same mission 
are often significant, if the statistical errors Oi are used to 
assess the significance. In what follows we compute the ad- 
ditional uncertainty introduced by the cosmic variance. 

If the true average XRB intensity is Ixrb then we want 
to compute the probability density function Pi{I\IxRB) of 
obtaining an intensity / in a measurement over a solid an- 
gle Qi with a statistical error Oi. This is computed by first 
convolving the confusion noise distribution with the statis- 
tical noise distribution (which we assume to be gaussian) 
and later shifting the resulting distribution to have a mean 
value Ixrb- This shift can be attributed to the behaviour 
of the source counts in the very faint regime, which would 
not affect the shape of the distribution. Assuming a top-hat 
beam (which is appropriate for our purposes) , the confusion 
noise distribution can be expressed as (Barcons 1992) 

Pi° nf (I) = J due' 2 ™ 1 exp(Q,i J dSN(S)[e 27Tt " s - 1])(1) 

where N(S) is the number of sources per unit solid angle and 
unit flux which in the range were we shall use it (S ~ 10~ 13 — 
1CP 12 ergcm~ 2 s _1 ) can be approximated as euclidean 

3 K ( S\- 5/2 

where we choose So = ICP 14 erg cm~ 2 s _1 and then K is 
the source density, extrapolated with the same euclidean 



law, at a flux limit 1CP 14 ergcm~ 2 s _1 . Indeed, the value 
of K depends on the energy band used for the observa- 
tions. For the 'hard' X-ray measurements [ASCA and Bep- 
poSAX) we assume an extrapolation of the Piccinotti et al 
(1982) source counts, now confirmed down to fluxes below 
1CT 1 ' 5 crgcm~ 2 s _1 (Cagnoni, Delia Ceca & Maccacaro 1998, 
Gendreau et al 1998). This gives K « 300 deg~ 2 . For the 
'soft' X-ray measurements we use the ROSAT normalisation 
K w 150 deg -2 (Hasinger et al 1993, Branduardi-Raymont 
ct al 1994). In modelling the confusion noise distribution 
for the various measurements of the XRB intensity, we have 
taken into account the possible excision of bright foreground 
sources by applying the corresponding cutoff in the flux 
integral of eqn. (1). This really only applies to the Chen 
et al (1996) ROSAT measurement where sources brighter 
than 3.4 x 10~ 14 erg cm~ 2 s -1 were excluded. For the re- 
maining measurements we have assumed that the sources 
in the Piccinotti et al (1982) catalogue had been avoided 
and therefore we applied a cutoff at a flux of 2-10 keV of 
3 x 10 -11 ergon -2 s _1 . 

Equation (1) ignores clustering of the XRB sources. In- 
deed clustering broadens the distribution of intensities with 
respect to a uniform distribution (Barcons 1992). From Bar- 
cons, Fabian & Carrera (1998) we see that for angular scales 
of the order of ~ 1 deg 2 , excess fluctuations due to clustering 
amount to ~ 0.5 — 1.2%, which implies a correction of the 
order of ~ 10% (at most) to the broadening computed by 
cq. (1). The relative amplitude of the clustering correction 
stays approximately constant when going to slightly smaller 
scales and decreases slowly when the small-scale declining 
part of the power spectrum of the fluctuations is reached. 
Given the smallness of this effect end the uncertainties in 
the modelling of the source clustering and its evolution wc 
decided to ignore this small effect. 

Convolving eqn. 1 with a gaussian of dispersion 
ai and shifting the mean to a value Ixrb gives 
the function Pi{I\IxRB) where now we express / and 
Ixrb as monochromatic intensities at 1 keV in units 
of keVcm -2 s _1 sr _1 keV -1 . For the conversion between 
monochromatic intensity and broad-band flux S per unit 
solid angle in a given energy band we use an XRB spectral 
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Figure 1. Bayesian density functions of the XRB intensity for 
the measurements obtained with ASCA, BeppoS AX and ROSAT. 
Individual measurements arc shown on the same curves as filled 
points (ASCA), asterisks (BeppoSAX) and squares (ROSAT). 
The HEAO-1 point is also shown as a vertical arrow. The top 
panel shows the bayesian distributions for each mission separately, 
the mid panel for the 3 missions all together and the bottom panel 
for the ASCA and BeppoSAX missions together. 



shape with a = 0.4 for the 'hard' instruments and a — 0.7 
for the 'soft' instruments which result in 

I(keVcnT 2 s" 1 sr" 1 keV" 1 ) = J- ^(deg 2 )" 1 A" 1 (3) 

where A w 200 if the flux is in the 2-10 keV band and 
A w 70 if the flux is in the 0.5-2 keV band. 



Typically, the curves including the cosmic variance are 
much broader than the the statistical error. To illustrate 
this, Table 1 also lists the dispersion E, obtained via a gaus- 
sian fit to Pi(I\IxRB) (which accounts for both the statis- 
tical and the cosmic variance) for comparison with the sta- 
tistical dispersion cr 4 . 

We now use these distributions to compute bayesian 
estimates of the true XRB intensity Ixrb- We assume an a 
priori distribution which is uniform over a sufficiently wide 
range of values of Ixrb (5 to 16 keVcm~ 2 s _1 sr _1 keV -1 ). 
The probability density function for Ixrb given a set of data 
is (Press 1989) 



niPi(ii\ixR B ) 

J dI XRB Hi Pi(Ii\I' XRB ) 



(4) 



where the products extend to the data points under consid- 
eration. In what follows, our Ixrb estimates correspond to 
the maximum of that function, which we also use to derive 
90% confidence intervals. 



3 THE INTENSITY OF THE XRB 

We first compute /(Ixrb) for the measurements obtained 
within each imaging X-ray mission separately. The results 
are shown in Fig. 1 where it can be seen that the various 
measurements with the same instrument are consistent when 
cosmic variance is taken into account. In particular the vari- 
ety of ASCA values which have been obtained with both the 
SIS and GIS instruments at different epochs, does not call for 
extra systematic effects. The same comment applies to the 
ROSAT data, points, as none of them is completely out of the 
distribution expected in terms of cosmic variance and statis- 
tical uncertainties, in spite of the fact that these intensities 
have been obtained with different PSPCs and/or different 
gains. The corresponding estimates of the XRB intensity 
are 9.8±?;5, ll.lt^ and 12.7^ keVcm" 2 s" 1 sr" 1 keV" 1 
(90 per cent errors) for ASCA, BeppoSAX and ROSAT re- 
spectively. 

Although there is considerable overlap between the dis- 
tributions corresponding to the three missions, there are ob- 
vious systematic differences. To illustrate this we combine 
all the data points with the use of eqn. 4 (see Fig. 1). The 
fact that there are points which obviously fall outside the 
distribution confirms the presence of systematic differences 
between the missions. Indeed the largest differences occur 
between ASCA and ROSAT amounting to ~ 30 per cent. 

Cross-calibration mismatches between missions have 
been reported by various authors, in particular between 
ROSAT and ASCA (e.g., Yaqoob et al 1994, Allen & Fabian 
1997). Iwasawa et al (1999) addressed this point with detail, 
motivated by simultaneous ROSAT and ASCA observations 
of the Seyfert galaxy NGC 5548. They find that the ROSAT 
PSPC spectrum is Aa ~ 0.4 steeper than the simultane- 
ous ASCA spectrum over a similar energy band. Iwasawa 
et al (1999) also comment on the fact that these calibra- 
tion mismatches do not occur between ASCA and any other 
missions, which find similar spectral shapes for a variety of 
objects observed. 

A further complication in the combination of XRB data 
from various missions is the different influence of the Galaxy. 
The main ASCA (Gendreau et al 1995) and BeppoSAX 
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(Vechhi et al 1999) measurements of the XRB intensity only 
used events above 1 keV, where the contribution from the 
Galaxy is ~ 1 per cent or less for these instruments. How- 
ever, in the Hasinger (1992) ROSAT estimate, the whole 
0.1-2.4 keV band was used to extract both the extragalac- 
tic and local components. For example, more than 30 per 
cent of the counts above 0.5 keV are of galactic origin. Al- 
though the influence of the Galaxy on the determination of 
Ixrb is uncertain in the ROSAT data, it is likely to be more 
important. 

Excluding the ROSAT data points, a bayesian estimate 
of Ixrb yields 10.0±g;| keVcm~ 2 s" 1 sr" 1 keV" 1 , but the 
fact that the largest measured intensity is the BeppoSAX 
one is suggestive of residual systematic effects (see again fig. 
1). In this case, however, these do not need to be larger than 
a few per cent. 

All of the above estimates give significantly higher 
values for the XRB intensity than the extrapolated value 
of the HEAO-1 data. For the BeppoSAX and the ASCA 
data, the HEAO-1 value (8keVcm~ 2 s _1 sr -1 keV -1 ) is still 
marginally consistent with the current modelling in terms of 
statistical and cosmic variance of the ASCA and BeppoSAX 
data. However, this interpretation would mean that the re- 
gions of the sky used in the imaging observations are system- 
atically brighter than the average sky sampled by HEAO-1. 
The possibility that the XRB spectrum between 1 and 3 
keV is steeper than the assumed a — 0.37 (there is a hint 
of this in the HEAO-1 data, see the mid panel of Fig. 4a in 
Marshall et al 1980) does not really help to alleviate this mis- 
match. The reason is that both ASCA and BeppoSAX mea- 
surements are quite sensitive to photon energies ~ 3keV, 
where the intensities found are significantly larger than the 
HEAO-1 result at the same energy. 



4 CONCLUSIONS 

Cosmic variance is able to account for the dispersion of the 
measured values of the extragalactic XRB intensity within 
the same mission. However, systematic differences remain 
among different missions, which cannot be understood in 
terms of spatial fluctuations. 

A cross-calibration mismatch between ROSAT observa- 
tions and those form other missions has been reported (Iwa- 
sawa et al 1999). But even using the ASCA and BeppoSAX 
data only the result (10.0±g;| keV cm~ 2 s _1 sr _1 keV -1 ) is 
higher than the XRB intensity at 1 keV extrapolated from 
the HEAO-1 data at photon energies above 3 keV. A steep- 
ening of the XRB spectrum at energies below 3 keV does 
not really help as both ASCA and BeppoSAX are sensitive 
to energies ~ 3keV and above where the mismatch persists. 
The only conclusion we can reach is that although internal 
calibration uncertainties in each mission amount to proba- 
bly less than 10 per cent, the cross-calibration among them 
has still large residual errors, preventing a more precise de- 
termination of the XRB intensity. 

XMM is and will be for many years the most sensitive 
X-ray imaging facility over the 0.1-10 keV energy band, pro- 
viding also moderate spectral resolution. Over 1 year, XMM 
will carry out observations of 500-1000 fields of ~ 0.2 deg 2 , 
so we expect ~ 100 deg 2 to be covered each year. The cos- 
mic variance (roughly scaling as Q -1 / 2 ) will then be small 



enough for unambiguous determinations of the XRB spec- 
trum and intensity. If the EPIC instruments can be cal- 
ibrated to significantly better than 10 per cent, this will 
certainly solve the issue of the intensity of the XRB. 
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